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By applying renormalization group method to the bosonized Hamiltonian of two-coupled chains 
with repulsive intrachain interaction, we have examined a role of backward scattering with a spin- 
anisotropy which competes with interchain hopping. From calculation of a dominant state in the 
hmit of low energy, it is found that superconducting state moves into spin density wave state when 
the anisotropy becomes larger than a critical value. Further phase diagram is shown on the plane 
of p-ology. 
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I. INTRODUCTION 

In low dimensional electron systems, which consist of 
an array of one-dimensional chains with repulsive inter- 
action, a noticeable phase diagram has been found where 
spin density wave |KDW) state competes with supercon- 
ducting (SC) statad. Among these materials, the organic 
conductors, (TMTSF)2^rjalts, have shown the SC state 
with d-wave like paringcrEl. Such a SC state is located 
close to the SDW state which exhibits one-dimensional 
fluctuation originating in intrachain interaction and a 
small interchain electron hoppingp. In addition to phe- 
nomena arising from the repulsive intrachain interaction, 
these materials exhibit anisotropy of spin susceptibility 
in the SDW state with the easy axis (^-direction) per- 
pendicular to the conducting chain although anisotropy 
is invisible above the SDW transition temperatureEl. 

For studying the competition between the SC state 
and the SDW state, we examine two-coupled chains with 
spin- anisotropic interaction. The electronic state which 
is obtained in a single chain with repulsive interactiora is 
quite different from those of quasi-one-dimensional sys- 
tem of these conductors in which the SC state shows 
d-wave like pairing denoting electrons between chains. 
Therefore, we need at least two chains treated correctly 
although infinite number of chains gives rise to a phase 
transition at finite temperatures. Actually such a SC 
state in case of weak coupling of Hubbard model has been 
obtained not only for two-coupled chains but also for 
three-coupled chains and speculated for arbitrary number 
of coupled-chaina]. It seems that the dominant state in 
two-coupled chains is relevant to electronic states for the 
above systems if we note that electronic states of quasi- 
one-dimensional conductors may be obtained as a result 
of developing these short range correlations perpendicu- 
lar to chain. 

Two-coupled chains with spin-isotropic interaction 
have been studied extensively. In case of forward scatter- 
ing described by Tomonaga-Luttinger model, the ground 
state given by the density wave is the same as that 



of a single chain except for rergaving degeneracy of in- 
phase and out-of-phase pairingsBilil. On the other hand, 
two-chains of Hubbard model with repulsive backward 
scattering leads tOj^he ground state given by the d- 
wave like SC statelBta'Ll, while that of a single chain 
is the SDW state. Numerical studies for the Hub- 
bard model also show an enhancement of SC state for 
some ranges of parameters of interaction, interchain hop- 
ping and fillingEsO. Analysis in terms of the bosoniza- 
tion indicates the fact that such a SC state appears 
due to formation of a gap of the total spin fluctuation, 
which is trisgered by the relevance of transverse charge 
fluctuationEHla. As for the spin-anisotropic interaction, 
the one-dimensional model has been exploped to study 
the SDW state of above organic conductorsta. 

In the present paper, the effect of spin-anisotropic 
backward scattering on the electronic state of two- 
couplecL-ehains is examined by improving the previous 
methodllS of the renormalization group, which is based 
on the bosonization. In section 2, formulation for renor- 
malization group equations is given where the backward 
scattering with parallel spins is distinguished from that 
with opposite spins to examine the spin-anisotropic in- 
teraction. In section 3, renormalization group equations 
are calculated numerically and the most dominant state 
in the limit of energy is examined where phase diagram 
of SC state and SDW state are shown on the plane of 
coupling constants and interchain hopping. Discussion is 
given in section 4. 



II. FORMULATION 

We consider the Hamiltonian for two-coupled chains 
given by 
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+9l^ 5p-p'5a-<j') 
^'^ki,p,aA^k2,-p.<T'.i^k2+q.-p',cr',i'^ki~q.,p',cr,i ■ (2.1) 

The quantity a\ ^ ^ ^ denotes a creation operator for 
the electron with momentum k, spin <j(=t, j or +,— ) 
where p = +(^) represents the right-going (left-g oing ) 
state and i{— 1,2) denotes chain index. In Eq. ( ^.1[ ), 
^k.p{= vp{pk — k-p)) is the kinetic energy with Fermi ve- 
locity vp and Fermi momentum kp. The energy, t, is that 
of the interchain hopping. Coupling constants gi\\ and 
gi± represent the intrachain interaction of the backward 
scattering with parallel and opposite spins respectively 
and g2 is the coupling constant of the forward scattering. 

By use of the unitary transformation given by 
Ck,p,a.ti = {-i^ak,p,c,,i + afe,p,cr,2)/\/2 (a* = ±), the first 
and the second terms in Eq. ( |2.lD are diagonalized 
resulting in the separation of the Fermi wave num- 
ber, fcp^ = kp — fit/vp. Based on the bosonization 
around the new Fermi point, we introduce the phase 
variables, 9^ and 0+ (0+ and 4>+), expressing fluc- 
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tuations of the total (transverse) charge density and 
spin density respectively. These variables are defined 

by E,^o2^Vj9^)e-"l''l/^^''?^4+,,p^.,,,Cfc,p,„,^ = [6+ + 

pd-+pi{0++p9^)+<j{<P++p(P^)+<jfi{(P++p(f>^)]/{2V2)E. 
They satisfy the relation, [9+{x),9-{x')] = 
= [9+{x),9_{x')] = [Mx)J-{x')] = 
i7rsgn(a; — x') where the suffix — denotes the canon- 
ically conjugate variable. The field operator is also 
expressed as '4>p.>j,p,{x) 
1/V27raexp [ipkp^x A- iQp 



,) exp {iirEp^^^i 



where 



e 



= [p9++9^+a{p(f,^ 



-n{p9++9-)+anip(l)^ 



0_)]/(2\/2) and a is of the order of the lattice constant. 
For keeping the anticommutation relation, the phase fac- 
tor, TrEp^a.ii, is chosen as S2„+j = A^i + • • • + A^2n + 
{-iy+\N2n+i+N2n+2)/2 with j^l,2 and n = 0, 1, 2, 3. 
The index i of Ni is defined hy i = {p, cr, /i) where 
(+,+,+) = 1, (+,-,+) = 2, (+,+,-) = 3, (+,-,-) = 
4, (-,+,+) = 5, (-,-,+) = 6, (-,+,-) = 7 and 
(—,—,—) = 8, respectively. Such a choice of Sp_cr,/i 
leads to a conservation of a sign of interactions. Thus 
Eq. (|2.l|) is rewritten as 
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27r2a2 



da; 



g„ cos — Atx/vpj cos'\/20+ + gb cos (V29+ - 4:tx/vpj cos\/20_ 

-f gc cos V26'_ cos \/20+ + gd cos V29- cos \/20_ 
-|- ge cos ^20+ cos |^\/26'-|_ — Atx/vp^ + 9f cos \/2(/'+ cos \f29- 

-I- 5g cos \/2(/)+ cos V20+ -I- gh cos V2(/)+ cos \/20_ 



(2.2) 



whe re vq = vp^ \ ~ {i^g2~ g^wM'^'^vvY , = 
«fv/1-{5iii/2^^^fP, Kg = [{l-(2g2~ffiii)/27rz;F}/{l + 
(252 - 5iIi)/2^«f}]'/', = [{1 + giii/27rz;F}/{l - 

5iii/27rwF}]^/^ and Kg ^ K^ ^ 1. In Eq. {^), coef- 
ficients are given by ga = .92 - Si 11, .9b = -^2, = 52, 
gd = -52 + ffiii and ge = -5/ = -gg = -gh ffi^ 
respectively. 

We examine the case where g^w = 52 but gi±^ is distin- 
guished from giw to examine the spin anisotropic back- 
ward scattering. By assuming the scaling invariance 
with respect to a ^ a' = ae'^', renorroalization group 
method is applied to response functionsE^lEll defined by 

^ G^(0 = -\ {Gl{l) My{l)) + G){1) + Gl{l) 
A G~{1) = \ {-Gl{l) My{l)) GUI) Myil)) 



Ra{xi -X2,Ti -T2) = {TrO a{xi,ti) 0\{x2,T2)j whcrc 

Tj is the imaginary time and Oa denotes operator for 
order parameter of SDW and SC states. In order to cal- 
culate second order renormalization group equations for 
all the coupling constants, we introduce g^, g~ and g~ 
given by 

K, = l + g„/2nvp {v^4>,9,4>) , (2.3) 

where = gi\\ and g^ — g-^ — 0. By making use of 
Gz — 5z/27n;F, renormalization group equations are ob- 
tained as 

cm , (2.4) 

Gl{l) + Gl{l)-gUl)My{l)) + gj{l)) , (2.5) 
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7) Jo(y(/)) + Glil) Jo(y(0) - G^(0 + gKI) - 6^(0 + gUI)) , 


(2.6) 








)Ga{l)~GS)GS) , 


(2.7) 




= 




)Gb{l)-G,{l)Gh{l) , 


(2.8) 




= {+G^{1) 


- 


)GS)-Gf{l)GS) , 


(2.9) 




= {+G^{1) 




)Gd{l)-Gf{l)Gh{l) , 


(2.10) 






- G~(0 


)Ge{l)~Ga{l)Gg{l)^Gb{l)Gh{l) , 


(2.11) 




= {-Gc^n: 


+ G~{1) 


)Gf{l)-G,{l)Gg{l)~Gd{l)G,,{l) , 


(2.12) 




= {-G4,{i: 




) Gg{l) - Ga{l) G,{1) Uy{l)) - GS) Gf{l) , 


(2.13) 




- {-GcbiK 




) Gh{l) - Gbil) GS) My{l)) - Gd{l) Gf{l) , 


(2.14) 


^^(^) 




{Gl{l) + 


Glil) + Glil)) Myil)) , 


(2.15) 



where the detail derivation is shown in Appendix A. In 
the above equations, i(/) — tH)/ ivpa^^), yil) = 4t(/) 
and Jnivil)) in = 0,1) is the n-th order Bessel func- 
tion. Initial conditions are given by 6*^(0) = gz/2'KVY 
and t(0) — t/ivpa~^). These equations retain the S U(2) 
symmetry for g^w — gi± — f/2- In Eqs. (2.7)-( 2.14 ), the 
bilinear term with respect to Gzil) iz — a ^ h) does exist 
for two-coupled chains with backward scattering. Such a 
term, whichJias been obtained by improving the previous 
calculationll3, is crucial to-reproduce the one-dimensional 
renormalization equationfl in the limit of small t. We ex- 
amine order parameters for the possible states in case of 
repulsive interaction, which are given by 



OlsdWi, 



Oss, 



a 

^aV^,.,^^'-,.,-^ , (2.16) 

a 

^(T/iV'+,a,A.'0- -<T,A. , (2.17) 



a.fi 



a.fi 
a 

= ^V+,a,MV'-,-a,M , (2.20) 

where ?/'p,<T,i(a;) = (l/VI) X^fc e'^'^afe^p^^^i. In Eqs. ( |2.16D - 
( pO| ), LSDW||,out (TSDW||,out) and CDWy.out denote 
longitudinal (transverse) SDW with intrachain and out- 
of-phase pairing and charge density wave (CDW) with 
intrachain and out-of-phase pairing while SSj^^in (TSj^_i„) 
represents the singlet (triplet) SC state with interchain 
and in-phase pairing. Since there is a symmetry with 
respect to gm = 2g2 in Eq. (^), the case with 2g2 > 
giw > is investigated explicitly where the relevant order 
parameters in terms of phase variables are reexpressed as 



TSDWi, 



Olsdw 



II .out 



OtsdW||. 
Oss, 



ix) 
ix) 
ix) 



'^exp -i{6l+ + cr0+}/\/2 cos |6'_ -I- cr0_ | /\/2 
'^fiexp -i{e++a4>^]/^/2 sin + cr(^+ 1 /\/2 
ex]i\i{e^+a(t)+]/V2 sin |L+ct^+|/%/2 



(2.21) 
(2.22) 
(2.23) 
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We study only 2A:f-SDW and SC states and do not con- 
sider 4fcF-CDW found for Kg < 1/2, of which denotes tl^ 
range beyond that of the conventional Hubbard modeC. 



III. SS STATE VS. SDW STATE 

We examine the state in the limit of low energy by 
calculating Eqs. (2.4)-( 2J^ ) where the corresponding 
energy is given by epe . The numerical calculation 
is performed with the fixed t/ep = 0.1 or 0.01 and 
several choices of coupling constants in the region of 
2^2 > > 0. The quantity eF(= vpa^^) denotes a 
Fermi energy and the coupling constant, which has the 
magnitude of th e o rd er o f the ban dwidth, is given by 
TTWp. Equations (2.4), (2.5) and (2.6) show that coupling 
constants corresponding to the total spin, the transverse 
charge and the transverse spin are renormalized leading 
to the excitation gap for the relevant case while that of 
the to tal ch a rge is unre norma lized . We find that, among 
Eqs. (|^, ( |2T2^ ) and (|^, LSDWJTSDW) state be- 
comes dominant for the relevant 0_ and </>_ (</>-, 0- 
and (f>+) and SC state becomes dominant for the relevant 
9_ and (h+. 




FIG. 1: The ^-dependence of Kj{l), K~{1), K^{1) and 
t{l)/e-p for t/ep = 0.1 and g]^ii/27rwF — .91^ /Sttwf — 
g^l^-Kvp = 0.4, where Ky{l) = 1 + Gu{l) {v = ^,0,0). 
The dotted curves denote K^il) of one- dimensional case 
given by i = 0. 

First the case with spin-isotropic interaction is shown 
in Fig. 1 where gi\\/2T:vp — gi±/2'!rvp = g2/2TTVp = 0.4. 
In this case, response function of LSDW is the same as 
that of TSDW for g^n = gi± due to the SU(2) symmetry 
for the spin rotation. The quantity t{l) is always rele- 
vant i.e. increases steeply for the present choice of pa- 
rameters. Renormalization of coupling constants shows 



that Gc{l) and Gf{l) and Gg{l) become relevant while 
Ga{l), Gh{l) and Ge{l) become irrelevant due to the spa- 
tially oscillating term in Eq. ( 2f2 ) and both \G d{l) / G c{l)\ 
and \Gh{l)/Gg{l)\ reduce to zero. Since K^{1) and K~{1) 

decrease and K~{1) increases, it turns out that the fluctu- 
ations for 0+, (/)_)_ and 0_ become relevant. The resultant 
excitation gap leads to a locking of phases at ^/2(l)+ — 0, 
=^0 and v^(^+ = TT (or ^20+ = tt, y/W- = tt 
and \/2(j)j^ = ) because the signs of the renormalized 
couphng constants for large / are given by Gc{l) > 0, 
Gf{l) < and Gg{l) > 0. Thus the dominant state in 
the limit of low energy is given by SC state with the 
singlet pairing of interchain and in-phase. Sue 
suit of SC state is consistent with the previous onell 
For parameters used in Fig. 1, SS (SDW) state becomes 
dominant compared to SDW (SS) state when the energy 
is smaller (larger) than a magnitude of ep exp[— 6.4]. For 
the comparison, the result for a single chain (i.e., t = 0) is 
shown by the dotted curve where Kj{l) = K~{1) — 1 and 

K^{1) decreases to 1. The difference between the solid 
curve and the dotted curve is attributable to the inter- 
chain hopping and appears for Z<^1, which corresponds 

to 4t(/)/eF>l. 




/ 

FIG. 2: The /-dependence of K~{1) with the fixed 
= 1, 1-02, 1.04 and 1.06 where t/ep = 0.1 and 
giii/27ruF = gil^-KVp = 0.4. 

Next the state for the spin-anisotropic interaction is ex- 
amined by varying the magnitude of the backward scat- 
tering with opposite spins, gi±. The case of gi± / gi\\ > 1 
is calculated where response function for LSDW becomes 
larger than that for TSDW. The case of 51^/5111 < 1 is 
discussed in section 4. In Fig. 2, the /-dependence of 
K~{1) is shown by choosing gi±/gii\ = 1, 1.02, 1.04 and 

1.06 where gi\\/2TTVp = g2/2'KVp = 0.4 and t/ep — 0.1. 
There are two kinds of fixed points depending on the 
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magnitude of gi^/giw- The quantity K~{1) with large 

I decreases for gi^/giw < 1.035 (case (I)) but increases 
otherwise (case (II)) although K~{1) with small I always 
decreases. In case (I), the behavior of renormalization 
is similar to that of the spin-isotropic case and then 
SC state becomes a dominant state. In case (II), the 
rapid increase of K~{1) shows the relevance </>_ instead 

of 0+ and then leads to the reduction of \G c{l) / G d{l)\ 
and \ G a(l ) I G h(l)\ which is opposite to the case (I). From 
Eq. ( 2.21 ), the dominant state in the limit of low energy 
for the case (II) is given by LSDW with intrachain and 
out-of-phase pairing since the signs of relevant couplings 
are given by Gd{l) < 0, Gf{l) < and Ghil) < and 
then phases are locked by -^720- 
V2^_ = (or 720- = TT, V2(/.+ 
In the onc-dimensional case,|_tl|. 
gi^ > is given by LSDWE3'E 

gap induced in the spin fluctuationE3. For two-coupled 
chains with gi± > g^w, the difference between response 
function of LSDW and that of TSDW is enhanced for 
I ^ I ln(4t*/eF)|, where tisjenormalized as t* due to one- 
dimensional fluctuationcJ^EJ. In the case (II), the domi- 
nant state is given by LSDW for all the energies while a 
crossover from LSDW to SC state with decreasing energy 
is obtained in the case (I). 



1.04 



= 0, V2(j)+ = and 
= TT and \/29- = n). 
dominant state with 
to the excitation 



1.02- = 




g/27IVp 

FIG. 3: Phase diagram on the plane of g and 5i-L/giii 
where where g = g2 = gi\\. The regions (I) and (II) 
correspond to SC state SDW state respectively and solid 
(dotted) curve denotes the boundary for t/e-p = 0.1 ( 
t/e-p — 0.01). In the inset, phase diagram on the plane 
of t and gi± / g^w with the fixed gjlixvp = 0.4 is shown. 

In Fig. 3, the ground state is shown on the plane of g 
and gi^/giii with the fixed t/e-p = 0.1 (solid curve) and 
0.01 (dashed curve) where g{= gm — 32). The regions 
(I) and (II) denote SS_L_in and LSDWy^outj respectively. 
For gi±/gi\\ = 1, we always obtain the SC state. For 
> Ij there is a critical value of g at which the 



SC state moves to the LSDW state with increasing g. 
Such a critical value decreases with increasing gi±/gii\ 
and/or decreasing t/cp- The inset shows a phase dia- 
gram of SC state (I) and LSDW state (II) on the plane 
of t and gi± / g^w with the fixed g/2TTVp — 0.4 where the 
SC (LSDW) state is obtained with increasing interchain 
hopping (anisotropy of interaction). For SC state, there 
is an energy gain by the ^c-term which competes with 
the ga-teim and then becomes relevant for energy smaller 
than t. On the other hand, the energy gain for LSDW 
comes from the g^-term, which is smaller but cooperates 
with the gd-term for all the energies. 



a. 
> 
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g2/27lVp 

FIG. 4: Phase diagram for two-coupled chains on the 
plane of 52 and g-^w where t/ep — 0.1 and gi± /gi\\ = 1.04. 
Regions (I), (II), (III) and (IV) denoteSS^ SDW, TS and 
CDW states which are given by Eqs. ( |2.17D , ( ^.16D , ( |2.20| ) 
and (2.1£) respectively. In the inset, the phase diagram 
of SS (I) and SDW (II) state on the plane of g and t in 
the case of 52 = ffiii =5 is shown where the solid curve 
and dashed curve denote the boundaries for gi± j g = 1.04 
and 1.02 respectively. 



In Fig. 4, the phase diagram is shown on the plane of 
g2 and g]\\ with the fixed g\^lg\\\ = 1.04 and t/ep = 0.1. 
The regions (I), (II), (III) and (IV) correspond to states 
for SS_L,in, LSD Wii.ou t, TS^n a nd C DW j_. out, w hich are 
given by Eqs. ( pT^ ), ( |2T^ ), ( |2^ ) and (pT9|) respec- 
tively. There is a symmetry given by gi — > g\\\ — 52 
with the fixed g-p^\ for the magnitude of response func- 
tions where response function of LSDW becomes equal 
to that of TS at g]\\ — 2g2. In the inset, the phase dia- 
gram of SS_L^in (I) and LSDW||^out (H) on the plane of g - 
t is shown for g-^i^ j g^w = 1.04 (solid curve) and g^±^ /g^ii = 
1.02 (dashed curve). Since these curves are well repro- 
duced by the formula, tjtp = exp[— Ci — C2 • 2'KVplg\^ 
with Ci ~ 1.3(1.2) and C2 ~ 2.6(1.2) for g^^jg^w = 
1.04(1.02), it is found that the boundary is determined 
by the competition between the interchain hopping and 
the spin gap induced by the spin-anisotropy interaction. 
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IV. DISCUSSION 

By applying the renormalization group method within 
the second order, we have investigated the dominant 
state in the hmit of low energy for two-coupled chains 
with the spin-anisotropic backward scattering. The 
phase diagram where SDW and SS states are calculated 
as the function of g^ii j 92 and t shows a noticeable inter- 
play of the interchain hopping and the intrachain inter- 
action. 

In the previous section, we have examined the case 
> 5iiij which leads to a competition between SC 
state and SDW state. For gi± < gi\\, we have found 
that results are always similar to those of Fig. 1 and 
then the dominant state is given by SC state. The 
difference between the case of gi± > gin and that of 
9i^ < originates in the fact that the renormaliza- 
tion for gi± > gm {gi± < gm) leads to a strong (weak 
coupling) coupling regime in a one-dimensional modeled. 
The case of gi± > gi\\, which leads to LSDW state, de- 
notes interactions which enhance the anisotropy along 
the direction parallel to the quantized axis. When the 
anisotropy is large for the perpendicular direction, one 
can expect TSDW in a similar way by adding interac- 
tions with the perpendicular anisotropy;, which have been 
treated in the one-dimensional modellS. |— . 

We compare the present result with that of FabrizioO 
who obtained the SC state for two-coupled chains of 
the Hubbard model by use of a diagrammatic method. 
The relation between our coupling constants and those 



of him is given by Gg 

■,(2) -X- „.(2) 



5" +5/" 
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(2) 



9f' 



^ 9' 
-9f^ 



- 9f' 

,(2) 



Gj = gll - .g/'l - 9^'^ 



(2) 



3/(2)^ G,^-2gJ, Gb 
2gJI, Ge = 2^6(1), Gf = -2ge 
2(7^^^ and Gh ~ ^Sg/*-^-*. He has derived eight 



Ga = 



- 9f 
.(1) 



equations by taking account of the SU(2) symmetry in 
the spin space, with the condition, = g^^^^ — gi^'^\ 
We have obtai ned eleven renormalization group equa- 
tions, i.e. Eqs. (2.4)-( 2.14 ), for coupling constants which 
become equal to his results within the second order un- 
der the assumption of the symmetry. Our renormaliza- 
tion group equations has been derived to cover the space 
where the spin-anisotropic interaction can be treated. 

Finally, we comment on the anisotropy of SDW state 
of organic conductors where the experimental evidence 
has been shown by spin susceptibilit}fl and the origin has 
been discussed in jtarms of spin-orbit coupling and dipole- 
dipole interactiontS. In the present paper, we have shown 
in Figs. 3 and 4 that a transition from SDW state to 
SC state occurs in the presence of spin-anisotropic back- 
ward scattering by the increase of the interchain hop- 
ping and/or the decrease of the renormalized coupling 
constant. Although such a transition might be expected 
under pressure resulting in the increase of band width 
and then vf, the critical value for {gi± — g^w )/ giw , which 
is evaluated by the moderate choices of parameters (i.e., 
g ~ ttwf and t/epi-p 0.1), is much larger than that found 
in the experimentQ. It is rather plausible to consider that 
the transitipn from SDW state to SC state in organic 
conductoraJ comes from the suppression of the nesting 
condition. However the results of the present calcula- 
tion indicate a reasonable magnitude of spin-anisotropic 
interaction for the competition between SDW state and 
SC state, which may be found in quasi-one-dimensional 
systems. 
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APPENDIX A: DERIVATION OF RENORMALIZATION GROUP EQUATIONS 



By treating the nonlinear terms in Eq. (2.2) as the per- up to the third order as 
turbation, the response function for field is calculated 
I 



^^^{t/V2),!,+ {xi.n)^-{i/V2)^+{x2,T2)^ 



-{K^/2)U{ri-r2) 



2 I d rgd r4 ^_(K^/2)U(r,-r2)^-2K^U(rs-r4) 



1 1 

'2! (47)2" J Q,2 Q,2 

X ^Q<^K^{U{ri~r3)-U{ri^r4)^U{r2-r3) + U{r2-ri}} ^ ^ 

X iG^e e cos2qQ(X3 — Xi) + Gje ' ' e' ' -f G^e ■> + Gj^e * 



1 1 

"3! (4^ 



,4 I ^Il^l±^Hp-(K4,/2)U{r,-r2)^-2K^U{r3-r^) 



e=±l 
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X exp 



GeGgC f> e * cos2qo(a;3 — X5) 

+6GbGeGhe e e f cos2(?o(a:3 - 2:5) 



'{2/K~)U{r3-r,,) -(2/ K~)U (rt-rs) ^ 



(Al) 



where U{r) = In 



I 

y^x'^ + VpT^ /a , (Pr = vpdxdT, x = the nonlinear terms, wiiJe these third terms are absent 
XI - X2 and r = n - T2. The quantity is replaced by one-dimensional caseB. These results originate in the 
VF and qo = 2t/vF. In order to obtain scaling equations nonlinear terms including both 0+ and 0± (or 0+ and 
of the coupling constants up to the second order, we need '/>± ) in Eq. (U) . By putting rs^r^ + r and = rs + r , 
the response functions expanded up to the third order of and expanding near r = 0, we find 
I 



Eq. (Al) = e 



-{K^/2)U{ri-r2) 



(47r) 



^\2 



g {K^/2)U{ri-r2)^-2K^U{r3-ri) 



X exp 



X G 



cff^ 



-2K~U{ra-ri) 



COS 2qo{x3 — X4) + G'j e 



-{2/K~)U{r3-ri) 



Gf~e 



-(2/K~)U(r3-ri)^ 



(A2) 



The quantities G° are given by 





— Gg 


— IGaGeGg 


Gf 


= G) 


— 2GcGfGg 


Gf 




— 2GaGeGg 




— 


— 2GhGeGh 



dr 



a \a 
dr 
a \a 
dr 



1-2K~ 



1-2K- 



1-2/ K~ 



/dr / r y 

dr /r\i-2/-ff~ 



— 2GdGfGh 



dr f r \ 
a \a/ 



1-2K~ 



a \a 

dr /r\i-2i<-' 



Jo(2gor)-2GcG/Gg 
'jo{2qor) - 2GrfG/G,, 



dr 



l-2/i<'~ 
e 

l-2/K~ 



(A3) 
(A4) 
(A5) 
(A6) 



r 



where r = {x^ + (wft)'^W^ and the second and third nentiating the third order terms of Eq. (Al). From Eq. 
terms of r.h.s. in Eqs. ([A^)-( A6) are obtained by expo- (|A2|), the quantity K^^ is derived as 
I 



Kf = K, '-GlKl 



- 2G)kI 



-.Gh^d 



/dr / r\ 
a \aJ 



Joi2qor)U{ri-r2) 



''U{r,-r2)--GlKl 



a \a 
dr / r \ 3-2X4,-2/ K- 



dr / r\3-2K^-2K~ 



(D 



U{r, - ra) 



C/(ri - ra) 



(A7) 



For a 



ae'^'El, these quantities are and G^*^(if^, G'^, gp, a') 



.2-Ki~K' 



scaled as Kf{Kl,g',q'^,a') = Kf{K^,g,qo,a), 
Gf{K,G',,q'o,a')^GfiK,,G,,q„,a) (a' /a) 
Gf{Kl, G'^.q'o, a') = Gf (if,, G„ go, «) (a'/a) 



.2-a:^-i/x~ 



Gf (if„G„go,a) 

(a' /a) * ^ * where iiT^, G^ and q'g denote renor- 
e , malized quantities. From this infinitesimal transform 



and putting K^{1) 



±1 _ 



1 ± G^{1), the renormalization 



2-K^~K~ equat io ns for the c oupli ng constants, Eqs. ( p. 41), ( p.ll[) 
Gg iK,G',,q'o,a') = G;; {K^,G,,qQ,a) (a' /a) * ( ^l2| ), ( |2l3| ) and ( pl|) are given by Eqs. 
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( A4), ( A5) and ( A6) respectively. In a similar way, renor- 
malization group equation for K2{1) {K~{1)) is calculated 



f rom the response function for 6 (</>) field, and Eqs. (2.7)- 
( 2.10 ) are obtained from the response function for both 

6 and (j) fields. 

The renormalization equation for t{l) is calculated as 
follows. The quantity. An, which denotes the actual dif- 
ference of the density between two bands, having kp^ 
and fcp- , is given by 



An = 2Aka 



T 



= 2{k-p+ ~ kY-)a + — j dx (It {k-p + ^ kp - j a , 

(A8) 

where {kp+ -kp_) / {2^ / L) x 2 = l/(2i) Ep,.,^ MPp,.,^. 
and Pp,a,^,{x,T) = ?Aj,,^^^(a;,r)'0p.a,M(a;,T). 

By using the phase variable, dxO+{x,T) = 
(7r/\/2)X]p,o-.pMPp,o-,^j(a;,r), the quantity An is ex- 
pressed as 

An = -2qoa + J dxdr (dxe+{x,T)'^ . (A9) 



By assuming the scaling relation for Eq. (A9) in a way 
simil ar to ref.tll, the renormalization equation for t{l), 
Eq. (2.15), is derived. 
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